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Geometrical constraints upon the unipolar model of V407 Vul and 
RXJ0806.3+1527 
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ABSTRACT 

V407 Vul and RX J0806. 3+1527 are X-ray emitting stars with X-ray light curves that are 
100% modulated on periods of 569 and 321 seconds respectively. These periods and no others 
are also seen at optical and infra-red wavelengths. These properties have led to the suggestion 
that the periods are the orbital periods of ultra-compact pairs of white dwarfs. There are two 
main double white dwarf models: the unipolar inductor model analogous to the Jupiter-Io 
system and the direct impact model analogous to Algol. In this paper we consider geometri- 
cal constraints on the unipolar inductor model, in particular what parameter values (compo- 
nent masses, orbital inclination and magnetic co-latitude) can describe the X-ray and optical 
light curves. We find that for a dipole field on the primary star, the unipolar inductor model 
fails to match the data on V407 Vul for any combination of parameters, and can only match 
RX J0806. 3+1527 if the sparser set of observations of this star have been unluckily timed. 
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1 INTRODUCTION 

In recent years much attention has been paid to what are possi- 
bly the shortest period binary stars known, V407 Vul (P = 570 
sec. lCropper et alll998tlMotch et alll996l) and RX 10806.3+1527 
(P = 321 sec. Ilsrael et alJll999h . Both of these stars show 
highly modulated X-ray light curves which are on for about 60% 
of the time and off for the remaining 40% I Croppe r et alJll998t 
Ilsrael et alj|l999h . Both stars also sho w only one period (and its 
harmo nics) at all wavelengths observed I R amsay et ali200oll2002t 
Ilsrael et all2002h . These and other properties have lead to the sug- 
gestion that the periods may be orbital, only possible for a pair of 
compact objects, most likely white dwarfs. This would make these 
systems strong emitters of gravitational waves and possible progen- 
itors of the semi-detached AM CVn stars. 

There are several rival models for these systems. V407 Vul 
was first suggested to be an intermediate polar (IP), in which case 
its 9.5 minute period would be the spin period of an accreting 
magnetic white dwarf <Motch et al]ll99fj) . However, the lack of 
any other period represe ntative of a longer period binary orbit 
lead ICropper et alj ll998h to suggest instead that V407 Vul was 



a polar containing a white dwarf with a strong enough magnetic 
field to lock its spin to the orbit, making it the first known "dou- 
ble degenerate" polar. This model received a blow wh en no po- 
larisat ion was detected iRamsav et alJl200oT) . which lead lWu et alJ 
(2002) to develop a unipolar inductor model. In the unipolar in- 
ductor model a slight asynchronism between the spin period of a 
magnetic white dwarf and the orbital period within a detached dou- 
ble white dwarf binary creates an electric current between the two 
components of the binary. The dissip ation of t his cu rrent powers 
the observed X-ray flux. Following IWu et alj |2002I>. an al terna- 
tive do uble white dwarf model was developed by Marsh & Steeahs 
1 2002). They proposed that accretion could be taki ng place without 
a disc iNelemans et aljEoOlklRamsav et alj|2002l) . since for peri- 
ods below ~ 10 min and for plausible system parameters it be- 
comes possible for the mass transfer stream to crash directly onto 
the accreting white dwarf. RX 10806.3+1527, although of a shorter 
period, is so similar to V407 Vul that all the above models apply 
equally well to it as they do to V407 Vul. Finally, the IP model was 
resuscitated by Norton et al. 1 2004) with a model in which we see 
the systems almost face-on with the X-ray variations caused by the 
accretion stream flipping completely from one magnetic pole to the 
other each cycle. 
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There is as yet no clear winner out of the various models, 
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all of which face difficulties. In the direct impact model, we ex- 
pect the (orbital) period to be increasing, whereas measurements 
show it to be decreasing in bo th V407 V ul l Strohmaverj20 l i 
2004 1) and RX J0806.3 +1527 lHakala et alj 120031: IStrohmav 
2003; Hakala et al. 2004). The weakness of optical emission lines 
in RX J0806.3+1527 and V407 Vul also seems hard to reconcile 
with an accreting binary, even though di rect impact m ay produce 
weaker line emission than disc accretion iMarsh & Steeehsl2002h . 
The weakness of the optical emission lines is probably also the 
most difficult fact to accommodate within polar or IP models since 
all such systems disco vered to d ate exhibit strong lines in their 
spectra. The IP model l Norton et al. 2004) is also unattractive for 
the fine-tuning it requires for two independent systems. The unipo- 
lar inductor model faces only a theoretical problem: it is short lived 
(~ 1000 years) because it lives off the spin energy of one of its 
two white dwarfs. As yet, no way of creating/maintaining asyn- 
chronism in the face of strong dissipation has been proposed. It 
does however nicely match the decreasing period of both systems, 
and so can justly b e claimed to be the leading model at present 
faakala et all2004 . 

Until a better alternative is found, we are faced with trying 
to select the "least worst" amongst the models, and every model 
must be tested as far as possible against the available observations. 
In this paper we examine whether the light curves are compatible 
wi th the unip olar inductor model, an area not considered in detail 
by Wu et al. 1 2002). I n section 3 w e define the geometry that we 
use to describe the Wu et al. 1 2002) model then we present our re- 
sults in section 4. We begin however by reviewing the observational 
constraints deduced from the X-ray light curves. 



2 OBSERVATIONAL CONSTRAINTS 

Our aim is to see if the observed light curves have a natural ex- 
planation under the electric star model. In doing so, it is important 
to avoid fine detail as probably all models fail at this level. For 
instance, in its simplest form of a small spot at the equator, the di- 
rect impact model gives 50:50 on/off, not the observed 60:40 ratio. 
However, it is not impossible to imagine that vertical structure near 
the impact site along with heating spreading downstream from it 
could lead to an increase of the "on" phase. Thus we try to dis- 
till key features which should be explained on any model, without 
being over specific in our selection. We identify the following key 
constraints: 

I. The bright phase during which X-rays are seen from the 
heated spots has to last for less than 0.6 of the cycle. While, as 
described above, it is possible to extend the visibility period of a 
given simple model, it is much harder to decrease it, so we believe 
that this is a fundamental restriction upon any model. 

II. The bright phase has to be more than 0.4 of the cycle because 
the maximum size of a spot (see later) extends the visibility period 
by 0.2 of a cycle at most. 

III. The phase of the X-ray maxima must deviate from a smooth 
long-term tre nd by at most 0.2 of a cycle, in order to match the 
observations ( Stro hmavej200ll2003h . 

IV. The light curves show no hint of eclipses and so we take the 
orbital inclination angle i to satisfy: 



Orbital axis 



cosz > 



ri + T2 




Secondary 

Figure 1. Orientation of the primary relative to the orbital axis and direction 
of the secondary centre. We also define the orientation of the magnetic field 
of the primary, tilted by fi from the z axis and whose projection in xy 
plane makes an angle tf> with the x-axis. The dot represents where the field 
line that goes through the secondary crosses the surface of the primary, the 
position of the spot makes an angle Q\ — 82 with the x-axis. 



3 THE GEOMETRIC MODEL 



(1) 



where r\ and are the radii of the two stars and a is the binary 
separation. 



Following the model presented bv lWu et alj 120021) we consider a 
system in which one star (which we call the secondary star) orbits 
within the magnetic field of its companion, the primary star. An 
asynchronism between the spin of the primary star and the orbit 
induces an electric current to flow along the field lines connecting 
the two stars, with energy dissipated mainly where these field lines 
enter the surface of the primary star. In the Jupiter-Io system the 
observed emission lies within 6° of the theoretically predicted foot 
point 1 Clar ke et alJll9961) . In our case we expect the emission to 
be even nearer to the foot point because of the smaller degree of 
asynchronism. 

Our task is to identify the positions of these foot points, and 
then work out for how long they are visible for a given orbital incli- 
nation. We treat them as points in order to allow the model to match 
the 40% dark phase during which X-rays disappear (constraint I, 
section |2j as easily as possible. To allow the model flexibility in 
matching the bright phase, we allow for the maximum size of the 
spots, as described at the end of this section. 

We will now describe the geometry of the model with refer- 
ence to Fig.Q We work with a right-handed Cartesian coordinate 
system (x,y,z) centred upon the primary star, with the a>axis point- 
ing towards the secondary star, the y-axis parallel to the direction 
of motion of the secondary star, and the 2-axis along the orbital 
axis (Fig. 0. The axis of the magnetic dipole is tilted by an an- 
gle /3 (often called the magnetic colatitude) with respect to the spin 
and orbital axes, which we take to be aligned. The azimuthal angle 
or magnetic longitude of the dipole ip is defined by the angle be- 
tween its projection onto th e orbital plane (the x-y-plane) and the 
s-axis. In lWu et alJ l2002l) 's model there is an asynchronism be- 
tween the period of the binary and the spin period of the magnetic 
white dwarf. Therefore ij) is the angle that completes a cycle every 
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few days owing to the asynchronism between the orbital frequency 
uj and the spin frequency of the primary star u s : 

ip — (lo s — LJ )t + ipo. (2) 

IWu et al ] 120021 estimate the asynchronism (ui s — u) )/uJo to 
be 1 part in 1000. This means that all possible azimuths of the mag- 
netic axis compared to the line of centres of the two stars are ex- 
plored within ~ 1000 orbits, which is only a few days. If the mag- 
netic axis is tilted with respect to the orbital axis, such azimuth vari- 
ation causes the phase of the light curve maxima to vary back and 
forth relative to any long term trend. By applying the III condition 
upon the phase, we effectively assume that the X-ray observations 
have sampled all possible relative orientations of the magnetic axis 
(all values of ip have been explored) and therefore we insist that the 
constraints of the previous section are satisfied for all ip. Ideally 
this would be the case, and indeed with sufficient observations one 
can get very close to it. We will consider how close we are to this 
in reality in section l4~Tl 

We define the positions of the foot points by where the mag- 
netic field line which passes through the centre of the secondary 
star crosses the surface of the primary star. This field line lies in a 
plane defined by the x-axis and the magnetic dipole vector which 
makes an angle a with the orbital plane Fig.Q We define the ori- 
entation of the foot points in this plane with two angles 0\ and 82 
which measure the angle between the dipole axis and a>axis and 
the angle between the dipole axis and the foot point respectively. 

Using the equation for dipole field lines r = C sin 8, we can 
deduce the following relation between 81 and 82 



(3) 



sint/2 = \ — sin e/i . 
a 



while from simple geometry one can show that 

cos 81 = cos ip sin j3, 



and that 



sin ip sin (3 



v/ (sin ip sin f3) 2 + cos 2 P 



(4) 



(5) 



To be able to calculate the visibility of the spot from the Earth, 
we need the orbital phase <p an d the orbital inclination angle i, 
which we define in the standard manner, i.e. i — 90° for an edge-on 
eclipsing system, = 0° when the magnetic primary star is at its 
furthest from Earth. With all angles defined, a given spot is visible 
while 



■*>)>-! 



where 



A = sin i y cos 2 (81 — 82 ) + cos 2 a sin 2 (8i —82), 
h = cos i sin a sin (#1 — 82), 



and 



tan 60 = 



cos a sin(#i — 82) 



(6) 

(7) 
(8) 

(9) 



cos(0i -0 2 ) 

For given values of f3, i, ip and n / a, it is a straightforward matter 
to calculate the remaining variables, and in particular to calculate 
the fraction of the cycle over which the spot can be seen. Before 
describing these calculations, we pause to consider the issue of the 
size of the foot points, which is important for the second of our 
constraints which defines the minimum period over which the spots 
are visible. 



The sizes of the foot points are related to the size of the sec- 
ondary st ar. If t he latter is tiny, then so too are the foot points. 
IWu et al] <2002l) give an equation for the size of the foot points 
(Al of their appendix A) from which we deduce that the maximum 
extent in longitude of the foot points, is of order 2r2/a radians. 
This will lengthen the bright phase by a corresponding amount. 
The maximum extent is thus set by the maximum relative size of 
the secondary, which is itself set by Roche geometry. For equal 
mass ratios, we deduce a maximum longitude extent of about 40°, 
which could extend the bright phase by ~ 10% of the cycle. To 
be conservative, we assume a maximum lengthening of 0.2 for all 
models, as outlined in constraint II of section|2| 

In our model we assume that the X-rays are only emitted at 
one of the two foot points. Allowing both spots to contribute makes 
matters worse. Although the restriction that the bright phase last 
for more than 40% of the cycle becomes easier to satisfy, it is more 
difficult to fit the 40% dark interval. The phase restriction is also 
affected negatively: for most configurations only one spot would 
be visible at a given time, but when the two spots are visible there 
would be a change in ip leading to a phase shift. 



3.1 Computations 

We computed the fraction of the cycle that one of the foot points 
is visible as a function of magnetic colatitude f3 (from to 180°) 
and orbital inclination (0 to 90°). In order to implement our as- 
sumption that we have observed all possible values of ip, we search 
for the maximum and minimum values of the spot visibility over a 
finely-spaced array of ip values from to 180° (symmetry allows 
us to avoid searching over 360°). In order to obey the constraints 
of section|2| we must have that 

(i) the maximum time that the spot is visible must be less than 
0.6 of a cycle, 

(ii) the minimum time that the spot is visible must be more than 
0.4 of a cycle. 

The masses of the white dwarfs have to be less than 1.4 Mq (Chan- 
drasekhar's limit) and larger than the Roche lobe filling mass: 
for RXJ0806.3+1527 Mi,M 2 > 0.12 M and for V407 Vul 
Mi, M2 > 0.065 Mq. We calculated the binary separation using 
Kepler's third law and the radius of the white dwarf using Eggle- 
ton's mass-radius relation quoted by Verbunt & Rappaport ( 1988). 



4 RESULTS 

We start by discussing how the constraints of section |2|restrict the 
possible values of the magnetic tilt (5 and the inclination angle i for 
typical white dwarf masses of 0.6 Mq (Fig.|2j- The first constraint 
(top panels) rules out combinations of low inclination and small 
angle between the magnetic field and the spin axis. This is because 
we would be looking at the north pole and the spot would also be 
in the northern hemisphere, i.e. it would be visible for most of the 
time which would violate the maximum 60% bright phase. Note 
that there is a critical magnetic colatitude when (3 C = 82 given by 



tan/3 c = y/ri/a 



(10) 



such that when the magnetic axis is tilted away from the secondary 
star (ip — 180°), the foot point lands upon the spin pole of the pri- 
mary star and will therefore always be visible, no matter the orbital 
inclination. For Mi = M 2 = 0.6 M , /3 C = 15° for V407 Vul and 
19° for RX J0806.3+1527. This critical value is clear in Fig.|5| 
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V407 Vul RX J0806.3+1527 




Magnetic tilt, |3 (deg) Magnetic tilt, |3 (deg) 

Figure 2. Geometrical constraints for the unipolar modal, for V407 Vul in the left hand panel and for RXJ0806. 3+1527 in the right hand panel. The shaded 
areas represent the parameter space ruled out by the restrictions I, II, III respectively in the top, middle and bottom panels. 
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Figure 4. The figure shows examples of X-ray phase light curves for 
V407 Vul. A, B and C represent different choice of parameters shown in 
Fig. with three possible -0=0°, 90° and 180° to illustrate the phase 
shift. For A, (3 = 5° and i = 87°, for B, /3 = 5° and i = 45° and for C, 
f3 = 135° and i = 87°. In all the figures the flux is normalised to unity. 



Figure 3. Track of the heated spot (dashed lines) viewed from above the 
primary star (outermost circle) for different choices of inclination of the 
magnetic field (3. For (3 smaller than the critical angle f3 c the track is always 
on the side facing the secondary (towards the right). In this case the X- 
ray pulse phase varies over a restricted range indicated here by the two 
tangential lines for f3 = 10°. For (3 greater than (3 C the heated spot makes a 
complete circle around the polar axis of the primary star causing a complete 
phase shift of the X-ray pulses over one beat period. 



The second constraint is symmetric to the first, so that the foot 
point cannot lie too close to the south pole without leading to too 
short a visibility fraction (central panels of Fig.|2j. 

The third constraint, that of phase, is highly restrictive as it 
removes all colatitude values in the range p c < f3 < 180 — j3 c . In 
other words, the magnetic and orbital axes must be nearly aligned 
in order for the phases not to wander. The reason for this is that 
for (3 in the excluded range, the phase switches by 180 degrees 
depending upon whether the magnetic dipole is oriented towards 
or away from the secondary star. This would show up as a large 
deviation in the X-ray pulses from a smooth trend. Note that our 



exact phase constraint (< 0.2 of a cycle variation) is actually a 
little more restrictive still since significant phase wander occurs as 
(3 nears (3 C . 

To better illustrate this, Fig. [5] shows the track of the heated 
spot on the primary star during one beat period as seen in the ro- 
tating frame of the binary. We can see that, for (3 < (3 C (or the 
symmetric case /3 > 180 — /3 C ), the phase offset is restricted to lie 
within the region defined by the two tangential lines to the small 
dashed circle but for f3 > {3 C the phase offset goes through a com- 
plete cycle. For the dipole field we assume, one can show that the 
tracks are perfect circles. 

To visualise the three constraints Fig.|4]shows examples of the 
light curves using some combinations of parameters that illustrate 
the problems of this geometry. When [3 = 5° and i — 87° (marked 
'A' on Fig. |2j the light curve obeys the first three constraints and 
it is similar to the observed data. In the second case considered 
(3 — 5° and i — 45° ('B' in Fig.|2j the spot is always visible and 
therefore fails to obey constraint I. In the third example /3 = 135° 
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and i = 87° ('C') the phase constraint is not complied with and 
there is a large phase shift between tj) = 0°, and ip = 180°. 

The fourth and last constraint is the restriction that the system 
does not eclipse. This is independent of magnetic tilt and depends 
only upon the masses of the two stars for a given orbital period. 

In figure [5] we show the possible values allowed by the com- 
bination of all the constraints mentioned above, for Ml = Mi = 
0.6 Mq, other white dwarf masses will be considered later. The fig- 
ures also show the maximum inclination for which the system will 
not have eclipses (horizontal lines). Together, the constraints rule 
out all parameter combinations. 

4.1 Detailed analysis of the phase constraint 

We assumed for simplicity that we have seen all possible values 
of the angle ip which measures the azimuth of the magnetic dipole 
within the binary. However, in this section we investigate further 
if this is really the case or if the observations of the two systems 
carried out so far are not sufficient to eliminate the possibility that 
the systems were only observed over a small range of ip, thereby 
increasing the available region of parameter space. 

We only used the X-ray data because according to the unipolar 
inductor model, the optical flux is from the secondary star and is 
there fore locked to th e orbit. 

IWu et alj (2002) deduced that for V407 Vul an asynchronism 
of 1 part in 1000 would be enough to explain the luminosity ob- 
served. 

The degree of asynchronism is a little arbitrary so we explored 
several values in the range 1 part in 100000 to 1 part in 700 for 
V407 Vul. We used the dates and durations of the X-ray obser- 
vations of V407 Vul to compute the phase shifts that would have 
been seen for a set of degrees of asynchronism in this range. The 
largest asynchronism of 1 part in 700 is set by the length of the 
longest observation and is such that this observation covered an en- 
tire beat cycle. A degree of asynchronism larger than this returns us 
to our original assumption of complete coverage of ip. The phase 
shifts were then fitted with a quadratic ephemeris and the scatter 
around the ephemeris was minimised by subtraction or addition of 
whole cycles where necessary. This process simulates the treatment 
of the observed pulse times by, for example. IStrohmavea (2004), 
and acts to reduce the phase shifts somewhat. Finally the values 
were minimised over the (unknown) phase offset ipo by repetition 
of the computation for 40 values equally spaced around a cycle and 
retention of the minimum shift calculated. This is in the spirit of 
trying to give the unipolar inductor model the benefit of the doubt 
where possible. 

We applied the same method to RX J0806. 3+1527. In this case 
the observations were much shorter and so our range of the de- 
gree of asynchronism is 1 part in 100000 to 1 part in 80 which 
corresponds to the period of the longest observation which was 
made with XMM-JVewton. We reduced archival XMM-Newton 
data from Novemb er 2002 to complement the phase timing by 
IStrohmaveil (200 3l) . Using the latest published value for the pe- 
riod P and its derivative P I Ha kala e t alj200 4|) and the p hase of 
the phase folded X-ray light curve of llsrael et ail fe003h MJD = 
52225.36153 ± 0.00004 for comparison, we obtained a phase shift 
for Nov 2002 close to zero, so we include the times for these obser- 
vations in our study. However, given present uncertainties in both 
P and P, we cannot rule out a larger phase shift between any of the 
observations. 

Our results are shown in Fig. [6j in which we plot the phase 
shift, minimised over ipo, as a function of the degree of asyn- 



chronism that would have been seen in the past observations of 
V407 Vul and RX J0 806.3+1527, for different inclinations of the 
magnetic field j3 if the lWu et alj i2002) model is correct. (Note that 
the phase shift is independent of the orbital inclination i.) 

As expected from our simple phase constraint, small angles of 
P such that (5 < j3 c do have low enough phase shifts to match the 
observations which show at most a deviation of 0.08 cycles for both 
systems. For V407 Vul, when j3 > Pc, even the minimised values 
that we plot are consistently larger than is observed, although the 
phase shifts do not match the value of 0.5 that we would predict for 
perfect observational coverage. In other words, the observations of 
V407 Vul are extensive enough that our conclusion that the region 
Pc < P < 180 — p c is ruled out by the absence of phase shifts 
remains unchanged. 

The much sparser observations of RX J0806. 3+1527 on the 
other hand, do not allow us to rule out large /3 angles with con- 
fidence as some frequencies show relatively small phase shifts. 
Clearly more observations of RX J0806. 3+1527 are required. Hav- 
ing said this, if the unipolar inductor model were to be let off the 
hook by this means, it would be at a significant price because then 
the X-ray phases measured so far would not necessarily be a true 
reflection of the orbital phase and the present indications of spin-up 
lHakala et alJi2004h . which depend in large part upon ROSAT data, 
could well be spurious. 



5 DISCUSSION 

In this section we will discuss the limitations of our analysis, the 
effect of varying the white dwarf masses and the effect of the spot 
size upon the allowed parameter space. 



5.1 Effects of different masses 

All of the constraints depend upon the masses of the white dwarfs, 
so in this section we explore whether changing the mass can allow 
some models to work. Increasing the size of the primary star helps 
by shifting the foot points toward its equator. Therefore, we expect 
an easing of the problem for small primary star masses, except that 
this will make the eclipse restriction harder to satisfy. Increasing 
the mass of the secondary star on the other hand is entirely bene- 
ficial since it increases the maximum inclination angle for which 
there is still no eclipse. However, it has to be said that there is some 
inconsistency here with respect to the foot point size which will be 
nowhere near our assumed maximum. If this could be taken into 
account reliably, it would count against large masses for the sec- 
ondary star, but we do not attempt to allow for this subtlety here. 

For V407 Vul we find that a very small area of the j3-i pa- 
rameter space is allowed for Mi = 1.4 Mq and M2 > 1.3 Mq, 
or for Mi = Mi = 1.3 Mq. The parameter space allowed re- 
quires the magnetic and spin axes to be nearly aligned (/? < 2° 
or P > 178°), the system to have a very specific orbital incli- 
nation (88.2 < i < 88.8°) and masses right at the top of the 
range for white dwarfs. Together the combination is very unlikely, 
and we will see below that it causes other problems. The case of 
RX J0806. 3+1527 is still tougher, and there are only possible so- 
lutions for Mi = 1.4 Mq and M2 = 1.4 Mq which are even 
more restrictive for both j3 and i (J3 < 1.5° or j3 > 178.5° and 
87.8 < i < 88.3° ) 

There are no known white dwarfs with masses as high as re- 
quired here. Moreover, even if the masses really were this high, 
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Figure 5. The lines going from top-left to bottom-right represent the parameter space ruled out by the combination of the constraints I, II and III. The only 
regions left out are two triangles in the top corners of the parameter space. However these are also ruled out when we include the absence of eclipses which 
rules out the regions shown by the oppositely sloped lines. Therefore we are left with no possible solution in the case of Mi = 0.6Mq and M 2 = O.GMq. 
Note that the vertical scale here has been adjusted to show only the top parts of Fig. 2 
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Figure 6. The minimum phase shift relative to a quadratic ephemeris as a function of the degree of spin/orbit asynchronism for /3 ; 
dashed line represents the maximum phase shift measured in the observations of the two stars = 0.08. 
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they would be in severe conflict with one of the main pieces of evi- 
dence in favour of the unipolar model, which are the measurements 
of decreasing period, presumed to be caused by gravitational radi- 
ation. We computed the values of the spin up of the systems for the 
lowest masses that have any available parameter space assuming a 
detached binary driven only by gravitational radiation loss. For the 
case of V407 Vul with Mi - M 2 = 1.3 M Q the spin up rate would 
be ~ 3.6 x 10~ 16 Hz sec -1 , so me fifty times the measured value 
reported bv IStrohmaved (2004) of 7.0 ± 0.8 x 10~ 18 Hz sec , 
while in the case of RX J0806.3+1527 for Mi = M 2 = 1.4 M Q 
the spin up rate would be ~ 2.3 x 10~ 15 Hz se c -1 , five times the 
measu red value reported by Strohmavei 1 2003) and Haka la et alj 
<2004h of 6.00 ± 0.1 x 10" 16 Hz sec" 1 . In consequence even the 
tiny region of parameter space that is opened up at high mass fails 
to match up against what we know of these systems. 



5.2 Spot Size 

We assumed a maximum spot size that would extend the visibility 
of a point spot for 0.2 of a cycle which corresponds to the sec- 
ondary star filling its Roche lobe , i.e. a very low mass secondary. 



But as seen in our analysis the problems facing the unipolar modal 
are eased at high masses. As the spot size depends greatly on the 
secondary mass, the actual size of the footprints for the masses con- 
sidered before would be only 0.02 of a cycle for Mi =0.6 Mq and 

0. 003 for Mi = 1.4 M in the case of V407 Vul. With such a small 
spot size it would be impossible to accommodate the model because 
for example if we force the bright phase to be more than 0.5 of a 
cycle, we eliminate all the values that were allowed by constraint I 

1. e. the first two constraints alone would rule out the model. 

There are some effects that might help make the spot larger, 
for example the X-rays could perhaps come from a vertical ex- 
tended region (along the field lines), that would allow us to see 
them for a little bit longer, although this might be difficult to rec- 
oncile with the thermal nature of the X-ray spectrum. Perhaps the 
foot points can leave a trail along the azimuthal direction (<f>o), that 
would enlarge the size of the emitting region. The theoretical esti- 
mates of the coolin g time scale indi cate that trailing is unlikely to 
be significant IStockman et alll994l) . but we can in any case show 
that increasing the visibility duration has only a limited effect in 
helping the unipolar model. In Fig.Qwe show the effect of increas- 
ing the spot size upon the available parameter space. Only for very 
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Figure 7. The top-left region of the parameter space where we show the 
constraint II for several maximum spot sizes 0.2 of a cycle (solid line), 0.3 
of a cycle (dotted line) 0.6 of a cycle (dashed line) for white dwarf masses 
of 0.6 Mq for 407 Vul. The allowed region is always above the line. So for 
the maximum spot size larger then 0.3 of a cycle there is a small allowed 
region in parameter space (triangle covering /3 > 173° 74° < i < 81°). 



large spot sizes, covering > 0.3 of a cycle in azimuth, is possible 
parameter space opened up. (Note that 0.6 is the largest possible 
azimuth extent given the need for a 40% off period.) 

In summary, even a fairly radical alteration to allow much 
larger spots than are predicted by the unipolar model at best opens 
up a small region of parameter space with (5 > 173°, 74° < i < 
81°. It is clearly unlikely that these applies to two independent sys- 
tems. 



6 CONCLUSION 

Assuming a dipolar field geometry and that the spin and orbital axis 
are aligned, we find that the unipolar inductor model presented by 
IWu et al] {2002) cannot match the X-ray light-curves of the two 
systems V407 Vul and RX J0806.3+1527. It fails because the foot 
points which produce the X-rays are situated quite close to the mag- 
netic poles, while at the same time the magnetic axes are forced to 
lie almost parallel to the spin axes of the primary stars to avoid ex- 
cessive phase shifts in the X-ray pulses. This then means that to ob- 
tain the near 50:50 on/off light curves observed, both systems have 
to be seen at such high orbital inclination (> 88°) that they would 
eclipse, and yet no eclipses are seen. High masses for the compo- 
nent stars permit tiny regions of viable orbital inclination/magnetic 
inclination parameter space, but lead to very large gravitational 
wave losses inconsistent with observed period changes. If for some 
reason the heated spots are much larger than the unipolar model 
predicts a small region of parameter space is allowed, but even this 
requires fine tuning for two systems. The only remaining chinks 
of light for the unipolar inductor model is that RX J0806.3+1527, 
owing to relatively sparse X-ray coverage, could have a relatively 
highly inclined dipole without our having spotted the large phase 
shifts so far, or possibly the field configurations are very different 
from a dipole in the two stars. Further observations are encouraged 
to close this loophole. We have been careful where possible to err 
in favour of the unipolar inductor model; we still have to conclude 
that it is not viable in its current form. 
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